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Abstract 

We calculate several differential distributions for exclusive double diffractive Xc(0 ++ ) production 
Q-f in proton-antiproton collisions at the Tevatron and in proton-proton collisions at RHIC and LHC 

jy^ ' in terms of unintegrated gluon distributions (UGDFs) within the fc^-factorisation approach. The 

uncertainties of the Khoze-Martin-Ryskin approach are discussed in detail. The g*g* — > Xc(0 ++ ) 
transition vertex is calculated as a function of gluon virtualities applying the standard pNRQCD 
technique. The off-shell effects are discussed and quantified. They lead to a reduction of the 
cross section by a factor 2-5, depending on the position in the phase space and UGDFs. Different 
models of UGDFs are used and the results are shown and discussed. The cross section for diffractive 
component depends strongly on UGDFs. We calculate also the differential distributions for the 
7*7* — ► Xc(0 ++ ) fusion mechanism. The integrated cross section for photon-photon fusion is much 
smaller than that of diffractive origin. The two components have very different dependence on 
momentum transfers t±,t2 in the nucleon lines as well as azimuthal-angle correlations between 
both outgoing nucleons. 
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I. INTRODUCTION 



The discovery of Higgs is the main motivation for the construction and putting into 
operation the Large Hadron Collider (LHC) at CERN. The analysis of inclusive cross section 
will be the "main road" of the future investigations. Different decay channels will be studied. 
The analysis in each (!) of these channels is rather complicated as huge unreduceable 
backgrounds are unavoidably present. 

The diffractive exclusive production of Higgs boson seems to be much cleaner in this 
respect. Many estimates of the corresponding cross section has been presented in the liter- 
ature. The so-called Durham model [l| is the state-of-art in this field. The cross section for 
the diffractive production is much smaller than the cross section for the inclusive case, but 
the ratio of the signal to more conventional background seems promising. Recently appeared 
a detailed analysis of diffractive production of MSSM Higgs [2J. We do not need to mention 
that any check of the used theoretical methods against experimental data are not possible 
at present, at least for the Higgs production. The way out is to study the diffractive pro- 
duction of heavy quarkonia where almost the same theoretical methods can be used. The 
basic diagram for the amplitude of the process is shown in Fig. HJ The production of heavy 
quarkonia received a lot of attention from both theory and experiment in recent years. For 
a review we refer to Refs. jl, 0, H, @, 0] ■ 




X h K lt HrjX 



FIG. 1: The sketch of the bare QCD mechanism. The kinematical variables are shown in addition. 

QCD dynamics at small quark and gluon momentum fractions (or large total energy), 
relevant for HERA, Tevatron, RHIC and LHC, is still poorly understood. It was shown in 
[§1. [To| that the combination of the fc^-factorisation approac h [III . [l~2l. 13 . 14] and the next- 



to-leading-logarithmic-approximation (NLLA) BFKL vertex [15] gives quite good agreement 
with data on inclusive QQ-P r °cluction. One can therefore hope that these concepts provide 
a valueable foundation also for exclusive processes. 

Kaidalov, Khoze, Martin and Ryskin proposed to calculate diffractive double elastic 1 
production of Higgs boson in terms of unintegrated gluon distributions [l|. It is not clear at 
present how reliable such calculations are, and it would be interesting to apply this approach 
to heavy quarkonia production. In Refs. [§, [l6| integrated cross section for exclusive double 
diffractive Xc(0 ++ ) production was estimated with the identical formulae as for the scalar 



1 Both protons survive the collision. 
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Higgs production with TH-* gg replaced by T Xc ( ++^ gg . Of course, such a procedure can 
be strictly right in general case only for ficticious structureless objects, when the internal 
wave function and gluon virtualities are neglected. From the spectroscopy point of view the 
X c (0 ++ ) meson is a quark-antiquark P-wave state, and it might be interesting to study an 
exclusive production of P-waves applying pNRQCD methods. Such a calculation may be 
especially important when we go to larger gluon virtualities. 

Parallel to the exclusive channel studies there was a lot of theoretical activity for in- 
clusive charmonium and bottomonium production; see for example [l~7l . 18 . 19| . There 



non-relativistic pQCD methods are usually applied. In these approaches the nonrelativistic 
quark-antiquark wave function is taken explicitly into the calculation. The vertex function 
g*g*qq corresponds to so-called quasi multi-Regge kinematics (QMRK), i.e. when q and q 
have similar rapidities and form a cluster. It is based on the formalism developed by Lipatov 
and Fadin [l5| for diffractive qq pair production. 

In the present work we shall use the pNRQCD methods for exclusive double diffractive 
X c (0 ++ ) production. The main ingredients of our calculations are the unintegrated gluon 
distributions (UGDFs) and the effective NLLA BFKL production vertex in QMRK (with 
reggeised gluon couplings to quarks). We would like to compare results obtained using this 
vertex with the KMR results. The projection of the heavy quark-antiquark pair onto the 
corresponding charmonium state is described in the standard way within the non-relativistic- 
quarkonium- model 0, 21, 22]. Finally, we shall refer to the KMR result [l6[ where the 



wave function is not included explicitly. For completeness, we shall include also photon- 
photon fusion mechanism of exclusive Xc(0 ++ ) production shown in Fig. [2J In addition to 
the diffractive QCD approach we discuss also a phenomenological Pomeron-Pomeron fusion 
approach. 




FIG. 2: The sketch of the photon-photon fusion mechanism. Form factors appearing in different 
vertices are shown explicitly. 

In the present paper we wish to calculate differential distributions for exclusive Xc(0 + ) 
production with different UGDFs from the literature. We shall calculate matrix elements 
for off-shell gluons. We shall discuss also uncertainty related to the approximation made, 
to the choice of the scale, etc. The contribution of the 7*7* fusion to the differential cross 
sections will be calculated. 
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II. DIFFRACTIVE QCD MECHANISM 



In the pNRQCD approach the diffractive exclusive reaction is viewed as shown in Fig. [31 
The corresponding calculation can be summarized as follows. First, the gg-continuum am- 
plitude is calculated. Then the gg — * qq amplitude is reduced to the gg — > Xc amplitude 
with standard projection techniques developed in [2(J 23 




FIG. 3: The basic diagram of double-diffractive charmonium production pp —* ppxc 



A. General kinematics of double diffractive QCD mechanism 

Let's consider first the gg-production. The kinematic variables for the process pp 
pp(qq) on the quark level is shown in Fig. HI 




FIG. 4: General kinematics of exclusive double-diffractive production. 

We use the following definition of the light cone coordinates 

k + = n+k a = k° + k 3 , k~ = n~k a = k° - k 3 , k t = (0, k\ k 2 , 0) = (0, k,0), 
where n are the light-cone basis vectors. In the c.m.s. frame 

P2 + Q0 - Pi- % 



n 



n 



(2.1) 



and the momenta of the scattering hadrons are given by 

vt=V%= y/s, Pi = p\ = Pi,t = P2,t = 0, 



where the Mandelstam variable s = 4E 2 ms . The momenta of the t-channel gluons are qo, qi 
and q2 (see Fig. H]). The on-shell quark and antiquark (with mass m) have momentum k\ 
respectively k 2 with 

, _ _ (™ 2 ~ k h) , _ _ (™ 2 ~ kit) 

In the high energy (large s) regime we have 

P + =qf -q£ « g+, P~ = gf - ?a ~ -^2", 

where P = k% + &2 is the momentum of the heavy charmonium with mass M: P 2 = M 2 ~ 
4m 2 . The longitudinal momentum fractions of the gluons are X\ = qf/pt, x 2 = —q 2 /p 2 - 
Of course, we will finally assume that ?o/i/2* = — 1*3.0/1/2,* | 2 - 

The decomposition of gluon momenta into longitudinal and transverse parts gives 

9i = xi(pi - q ) + qi, t , Q2 = ~x 2 (P2 + Qo) + q2,t (2.2) 
< £1,2 < 1, q = x'iPi + qo,t = -x' 2 p 2 + qo,t • 

We take into account below that x[ = x' 2 = Xq. Making use of conservation laws we get 

Qi + Pi =Pi~ Qo, Q2+P2 + qo= P 2 ■ (2.3) 
Taking the transverse parts from these relations gives 

qi,t = -{p[,t + qo,t), q2,t = P 2 ,t - Qo,t ■ (2.4) 



B. Matrix element for exclusive double diffractive x c (0 ++ ) production 

According to Khoze-Martin-Ryskin approach (KMR) we write the amplitude of ex- 
clusive double diffractive colour singlet production pp — > ppXcj as 



M - 2^ 2NT^1^ ] * % > tVj qj^qf, ^ 

The normalization of this amplitude differs from the KMR one [l|, [l6[ by the factor _sj 2 and 
coincides with the normalization in our previous work on exclusive ^'-production [24j. The 
amplitude is averaged over the colour indices and over two transverse polarisations of the 
incoming gluons [l|. The bare amplitude above is subjected to absorption corrections which 
depend on collision energy. We shall discuss this issue shortly when presenting our results. 

The vertex factor Vj lC2 = Vj lC2 { y q\ v q\ v Pl /l ^) in expression (I2.5P describes the coupling 
of two virtual gluons to Xcj- m eson that can be written as 

Vj ClC2 = V(qq -> Xcj) • n° 2 (ki, h) , (2.6) 

where V(qq — > Xcj) is the operator that projects the qq pair onto the charmonium bound 
state (see below), ^ ClC2 (ki, k 2 ) is the production amplitude of a pair of massive quark q and 
antiquark q with momenta hi, k?, respectively. 
Within the QMRK approach [l5| we have 

*( Cl , c 2 ; i, k; k u k 2 ) = - g 2 ( t %t? k b(k u k 2 ) - t c k )tp(k 2 , k,)), a s = £ , (2.7) 
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FIG. 5: The effective vertex within QMRK 



where t c are the colour group generators in the fundamental representation, b, b are the 
effective vertices (12.81) arising from the Feynman rules of QMRK illustrated in Fig. (jSJ). 



Kfa, k 2 ) = 7 
b(ki, k 2 ] 



qi - h 



m 



7 



(qi — kx) 2 — m 2 

- qi - h + m 
(qi — ki) 2 — m 2 



T- 
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7/3 r+ (g2,gl) 

{h + k 2 ) 2 

7 /3 r + - /3 (g 2 ,gi; 

(h + k 2 ) 2 



(2.8) 



Taking into account definitions (12.1 1) and (12. 3p and using the gauge invariance property 
we get the following projection 



ycxci 



<n- v V™{ qi ,q 2 ) = -^ q -Mv^( qi ,q 2 ). 



(2.9) 



Normalization of the polarization vectors coincides with the one in Ref . [25J . Since we adopt 
here the definition of the polarization vectors proportional to gluon transverse momenta 
qi/2,ti then we must take into account the longitudinal momenta in numerators of vertices 
( 12. 8p . While projecting on the color singlet the ggg-verticies T H in (12.81) cancel each other 
and disappear from the resulting matrix element, so there contributes only the first diagram 
in the decomposition in Fig. [51 

As a consequence of the gauge invariance the charmonium production amplitude for on- 
mass-shell quark and antiquark states obeys the following property 

u{k 1 )^! c ^v{k 2 ) -> if q u or q 2>t -> 0. 

Projection of the hard amplitude onto the singlet charmonium bound state V^ C2 is given 
by a 4-dimentional integral over relative momentum of quark and antiquark q — (k± — k 2 )/2 



Vj^(quq2)=V(qq^ Xc j) 



{fa, h 



2tt 



^ *—t Jm 



i,k L Z ,S Z 



d 4 q J o q 2 



x 



X 



$ w ,(q) • (L = 1, L z ; 3 = 1,S Z \J, J z ) (3z, 3fc|l) Tr {^TJPs^,} 



(2.10) 



tyClC 2 
IK, [IV 



-9 



J.C1J.C2 

Hj L jk 



qi ~ fa 



m 



(Qi — ki) 2 — m 2 



7u 



_ J.C2 J-Cl 



7*. 



Qi - k 2 + m 

(<7l - ^2) 2 - 
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Here the function $>l=i,l z {q) is the momentum space wave function of the charmonium, and 
the projection operator Vs=i,s z f° r a small relative momentum q has the form 

Vs=i,s z = 7^(k - m)^-(k + m) . (2.11) 

Factor 2ir in the formula (12.101) has been introduced to compensate the difference be- 
tween the analogous formula in Ho, Q written in 3-dimentional form. The Clebsch-Gordan 
coefficient in color space in our case is (3i,3k\l) = 5 th / ' \fW c , where factor l/\/W~ c provides 
the averaging of the production matrix element squared over intermediate quark states in 
the loop. Using that we get 

rp/ lTf ^ v XCl c 2 eP ( S z) 9 2 m f f qi-k-m gi - k 2 + m \ 

Tr(*V s=1 , Sz ) = -5 7^^Tr| ^- - ^ _ 7 „ - 7 „_ - ^ - _ 7 „j x 

x (fc 2 -m) 7p (fci + m)j. (2.12) 

Since P-wave function $l=i i l z vanishes at the origin, we may expand the trace in (12.101) in 
a Taylor series around q = 0, and only the linear terms in q a in the trace f 1 2 . X 2 j) survive. 
This yields an expression proportional to 



/ |^k $ ^ (q) = -Vi eCT(L ^ m (2 - i3) 



with the derivative of the P-wave radial wave function at the origin 1Z'(0) whose numerical 
values can be found in [26]. The general P-wave result ( 12.101) may be further reduced by 
employing the Clebsch-Gordan identity which for the scalar charmonium Xcj=o reads 



T/4 ee £<1, L z ; 1, S,\0, 0) e°(L z )e»(S z ) = ^| 



ptrpp 

M 2 



Using the relations (I2.3P and q\ — q 2 = Pm an d s 3> |qo,t| 2 , we obtain 

= \{qlt + Qtt - M 2 ), sx lX2 = M 2 + |P A d 2 = M| , (2.14) 
that will be useful below. The longitudinal momentum fractions are now calculated as 

v/M 2 + |P Mit | 2 

xi i2 = = exp(±y) , 

V* 

x 'i,2 = x o = -^7^- • (2.15) 

Above y is the rapidity of the produced meson. 

Therefore, within QMRK approach we get finally the following vertex function 

Vj=o{quQ2) = 8ig 2 



M ^MW C (M'-qh-qh) 2 



7 



We have also calculated the subprocess matrix element squared B(qi,q 2 ) = V IW V 11V that 
is usually used in inclusive production calculations. The form of B(qi, #2) is identical to the 
form of the matrix element squared obtained very recently by Likhoded and Luchinsky in 

Ref. Ur 



The objects f° g s J (x u x\, gg >t , q\ p t x ) and f° f 2 f (x 2 , x' 2 , gg jt , g| t , t 2 ) appearing in formula flZSJ) 
are skewed (or off-diagonal) unintegrated gluon distributions. They are non-diagonal both 
in x and q\ space. Usual off-diagonal gluon distributions are non-diagonal only in x. In 
the limit xx j2 — > x' 12 , ql t — > q\i 2t and ti }2 — > they become usual UGDFs. Our choice of 
different UGDFs will be discussed in more detail in a separate section. 



C. Khoze-Martin-Ryskin approach 

In the original Khoze-Martin-Ryskin (KMR) approach [l[ the amplitude is written as 
M = N f d2qo f i X f +)] f™*{ Xl , x[, Ql, ^ h)ff MR {x 2l x' 2 , Ql„ t 2 ) , (2.16) 

J %,t1l,tQ.2,t 

where only one transverse momentum is taken into account somewhat arbitrarily as 

Ql,t = mm {9o, f > it} > Ql,t = min {9o,i> it) , ( 2 - 17 ) 

and the normalization factor iV can be written in terms of the Xc(0 + ) — > gg decay width 
(see below). 

In the KMR approach the large meson mass approximation M ^> |qi,t|, |q2,t| is adopted, 
so the gluon virtualities are neglected in the vertex factor 

P[ X c(0 + )] ~ {q u q 2tt ) = (q ,t + p[, t )(qo,t ~ P 2 , t )- (2-18) 

In our approach we generalize the approximation taking into account corresponding off-shell 
effects. 

The KMR UGDFs are written in the factorized form: 

f« MR (x,x',Ql^;t) = f g KMR (x,x',Ql^)eMbot) (2.19) 



with bo = 2 GeV~ 2 [1]. In our approach we use different parametrization of the t-dependent 
isoscalar form factors (see Eqns. (I2.32j) and (I2.33|) below). 

Please note that the KMR and our skewed UGDFs have different number of arguments. 
In the KMR approach there is only one effective gluon transverse momentum (see Eq. (12. 171) ) 
compared to two idependent transverse momenta in our case (see Eq. fl2.32p ). 

The KMR skewed distributions are given in terms of conventional integrated densities g 
and the so-called Sudakov form factor T as follows: 



f« MR (x,x',Qlv 2 )=R : () 



T(Ql^)xg(x,Q 2 t ) 



[2.20) 



The square root here was taken using arguments that only survival probability for hard 
gluons is relevant. It is not so-obvious if this approximation is reliable for cc quarkonium 
production. In addition this has to be contrasted with the choice of gluon momentum of the 
KMR UGDF in (12.171) as minimal (not harder) of two gluons. The factor R g approximately 
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accounts for the single logQ 2 skewed effect Please note also that in contrast to our 
approach the skewed KMR UGDF does not explicitly depend on x' (assuming i'<Ci< 1). 
Usually this factor is estimated to be 1.3-1.5. In our evaluations here we take it to be equal 
1 to avoid further uncertainties. 

In contrast to the Higgs case, in the case of light quarkonium production rather small 
values of gluon transverse momenta give the dominant contribution to the integral (12.161) . 
Therefore it becomes essential what one does phenomenologically in the nonperturbative 
region Q 2 < Qq, where the Ql is a minimal nonperturbative scale for standard integrated 
distributions. Of course, formally for scales smaller than Qq the standard collinear gluon 
distributions do not exist and an extra extension is unavoidable. This issue was not discussed 
in detail in the literature. We shall illustrate this point in the result section. 

The Sudakov factor is the result of resumming the virtual contributions in the DGLAP 
evolution and reads 



T{Ql " 2 



// 



exp 



a s {k 2 



2tt 



dk 2 



k 2 



1-A 



(2.21) 



with A = k t /(fi + k t ). In their (KMR) estimates the hard scale is usually taken // = M x /2. 
Of course the choice of the scale is somewhat arbitrary, and the consequences of this choice 
were not discussed in the literature. 

In Figj6] we show the KMR distribution (12.201) as a function of Qf for different values of 
x specified in the figure for two different choices of the scale. The DGLAP distribution for 



Qt < Qo * s n °t wen defined and we arbitrarily put it to zero as in Ref.[16jj. In principle, one 
could try several other extrapolations into the nonperturbative region to see its influence on 
the resulting differential cross sections. 
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6: KMR distribution as a function of effective transferse momentum Q\ for different values 



M 2 (left panel) and fi 2 = (M x /2) 2 (right panel). 
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D. Our vertex versus Khoze-Martin-Ryskin vertex 

We wish to compare our vertex function (12.161) with the KMR vertex function [l6|. In 
the KMR limit of large meson mass M ^> |qit|, |q2,t| and SX1X2 — M 2 the vertex reads 

Following now the KMR notations Q we write the total NRQCD QMRK amplitude (1231) 
(averaged over colour and polarisation states of incoming gluons) in the limit M 3> qij, Q2,t 
as 



A1 = Att - / d go,t-P[Xc(0 + )J - 2 — 2 — 2 ' ( 2 - 23 ) 



where the normalization is 



and the vertex factor is defined in (12. 18[) . 

Normalization constant A can be obtained in another way in terms of the partial decay 
width r(x c o — > gg) as (see formula (19) in the KMR paper [l6l |) 

= K (Nf^M* ' (2 ' 25) 



Using the expression for r(x c0 — > gg), obtained in the framework of pNRQCD [27], 128 



I7£'(0)l 2 

r(Xco - gg) = 32N c a 2 J-^^, N c = 3 (2.26) 



we get the normalization constant (with K — 1) 



A=ig i.,*m (2 . 27) 

which coincides with the normalisation of the vertex factor obtained within the QMRK 
approach (I2T241) . 

Therefore, in the leading order the QMRK approach is in agreement with the KMR 
approach in the limit of large meson mass M ^> |qi,t|, |q2,t|- We shall discuss deviations 
from this approximation due to the off-shell effects in the result section. Similar analysis of 
off-shell effects was performed recently for inclusive Higgs production in Ref. (29| . 



E. Off-diagonal UGDFs and choice of scales 

In the present work we shall use a few sets of unintegrated gluon distributions which aim 
at the description of phenomena where small gluon transverse momenta are involved. Some 
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details concerning the distributions can be found in Ref. |30J. We shall follow the notation 
there. 

The larger energies, the smaller values of parton momentum fractions come into game. 
Therefore at larger energies we shall use distributions constructed exclusively for small values 
of x. Two of them are based on the idea of gluon saturation. One of them was obtained based 
on a saturation-inspired parametrization of the dipole-nucleon cross section which leads to 
a good description of the HERA data [31]. The second one 32| was constructed to describe 
the inclusive RHIC pion spectra. The third one is the asymptotic BFKL distribution 33 . 
We do not wish to repeat more details here. It can be found in individual references as well 
as in Ref. 3(| where applications of UGDFs to cc correlations was discussed. 



Due to its simplicity the Gaussian smearing of initial transverse momenta is a good 
reference point for other approaches. It allows to study phenomenologically the role of 
transverse momenta in several high-energy processes. We define simple unintegrated gluon 
distribution 

^ auss {x, kl /4) = xg coll (x, /4) • f Gauss {k 2 t ) , (2.28) 

where g co11 (x , fi 2 F ) are standard collinear (integrated) gluon distribution and fcaussik 2 ) is a 
Gaussian two-dimensional function 

f G auss(k 2 ) = -L exp {-k 2 /2a 2 ) /vr . (2.29) 

The UGDF defined by Eq. fl2T28|) and (12T29|) are normalized such that 

J ^ auss (x, k 2 , fi%) dk 2 = xg cM (x } /4) . (2.30) 

The UGDFs have the following property 

f(x,k 2 )^0, (2.31) 

if A;^ — ^ 0. The small- k 2 region is of nonperturbative nature and is rather modelled than 
derived from pQCD. 

The two-scale off-diagonal distributions require a separate discussion. In the general case 
we do not know UGDFs very well. It seems reasonable, at least in the first approximation, 
to take in the amplitude (12.51) 



/;f = \/n i '(x' 1 , q i t ,^)-n i> (x 1 , q i t ,^)-F 1 (t 



& = ^ff\A,qlt^l)-ff\^ql t ^ 2 ) ■ F 1 (t 2 ) , (2.32) 
where -Fi(ti) and Fi(t 2 ) are isoscalar nucleon form factors [3~ij 



Fi{ti,2) — 7-; — o — 1 — ui 1 — 7nvT\2 ' (2.33) 

(4m^ - ti, 2 )(l - ti, 2 /071) 2 

and t\ and t 2 are total four-momentum transfers in the first and second proton line, re- 
spectively. The proton form factor in the form (12.331) gives rather good description of the 
t-dependence of the elastic pp cross section at high energies, i.e. for kinematics similar as in 
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our case. The above prescription for UGDFs (I2.32p is a bit arbitrary, although it is inspired 
by the positivity constraints for collinear Generalized Parton Distributions [35j. It provides, 
however, an interpolation between different x and qf values. Our prescription is more sym- 



metric in variables of the first and second exchange than the one used in Ref. |36J for Higgs 
boson production. 

The choice of the (factorisation) scales here is not completely obvious too. We shall try 
the following three choices: 

(1) ^ = M 2 , f = M\ 

(2) ti = Ql /i 2 = M 2 , 

(3) fi 2 = g 2 1 (+freezing at g 2 1 < Q 2 ), /i 2 = M 2 . (2.34) 

The first choice is similar as in Ref. [l|, 37]. However, it is not obvious if the scale 
associated with the "hard" production (g*g* — > Xc) can be used for the left part of the 
gluonic ladder where no obvious hard scale appears. Therefore we shall try also the second 
choice where we shall use Q$ = 0.26 GeV 2 , i.e. the nonperturbative input for the QCD 



evolution in Ref. [38] . Another option was proposed by Lonnblad and Sjodahl in Ref. 
They take g 2 1 as a first scale. In our case this prescription must be supplemented by freezing 
the scale for gluon transverse momenta smaller than Q (minimal perturbative scale). 

When inspecting Eqs. ( 12.51) and (12.321) it becomes clear that the cross section for elastic 
double-diffractive production of a meson is much more sensitive to the choice of UGDFs 
than the inclusive cross sections. 



III. 7*7* FUSION MECHANISM 

As stated in the introduction we wish to investigate the competition of the diffractive 
QCD mechanism discussed in the previous sections and the 7*7*-fusion mechanism shown 
in more detail in Fig. [7J 




FIG. 7: Kinematics of exclusive 7*7* fusion mechanism of Xc-meson production. 
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A. NRQCD approach 

In the most general case the Born amplitude reads 



M XlM , K , x , = ±e^u(p[,\[)[F 1 (t 1 )Y± i ^qiy'F 2 (t 1 )}u(p 1 ,\i)} ^ Vff^ J (h,k 2 ) ^ 



9w 



a 



{u(p' 2 , X' 2 )[F x {t2)^ ± i^-qi,^F 2 (t 2 )]u(p 2 , A 2 )} , 



(3.1) 



where the sign " + " stands for pp-scattering and the sign " — " for pp-scattering. In the 
following we have omitted the spin-flipping contributions related to the respective Pauli 
form factors. Limiting to large energies (y/s ^> M + 2m7v) and small momenta transfer t\ y2 
(|*i,2| ^ 4m^r) the matrix element for pp — > ppXcj reaction via 7*7*-fusion can be written 

as 



M 



77 



eF l{ t x )^±^V::^^ q2 ) ^±^eF x ii 2 ) 



(3.2) 



where i'i(ti) and Fi(t 2 ) are Dirac proton electromagnetic form factors, and the 7*7* — > Xcj 
vertex has analogous form as f)2.10p 



yi 7 -»Xcj( 



/1, fc) = V(qq - Xcj) • <M*i> y= 2 ^EE^/ (?° - |) >< 



L Z ,S Z 



X 



$ L=1 , Lz (q) • (L = 1, L z ; 5" = 1, ^| J, Jg) (3z, 3*|1) Tr {*%V s= i, St } 



(3.3) 



. Op 

v I 3 



q\ — h\ — m 
(qi — ki) 2 — m 2 



IfJ, In 



qi — k 2 + m 
(qi - k 2 ) 2 - m 2 



(3i,3A;|l) 



and 6 Sik = N c . Repeating similar steps as for diffractive production we get 



yi 7 -*Xcj{ 



2(?2,v - q\,u)q\, P 9na - 



{qi,h - q2,p,)(q2,a + qi,a)gv P + -(qi, P - <? 2 , P )(<?i, CT - <h, a )g^ + (M 2 + 3? - q\)gv P g^ + 



+ 



qi,a + <?2,<x 



{qi, P - q2, P ){qi,^qi,u - QVQV) + 



2(gig 2 ) 2 

+ (M 2 - q\ + ql)qi, p g^ ~ 2M 2 q 1 ^g up + {q u fi «-> g 2 , ^} 



(3.4) 



For example, in the case of scalar meson (J = 0) we have 
Vp^ Xc0 (qi,q2)- 



. 2e\ 2 ft'(0) M 2 g lxv - (g M - g 2 , M )(g v - g 2 ,„) 

1 3 J M VvrM\ 



[M 2 (M 2 - 2g 2 - 2g 2 ) + {qi - g 2 ) 2 ]^ + 2(M 2 + qf - g 2 )g 2i/1 g v + 



%ig 2 ) 2 . 

+ 2(M 2 + q\ - qf)qi, u qi,n + ^M 2 q 2 , v q^ 



(3.5) 
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This amplitude can be rewritten in the general gauge invariant form 2^] in terms of two 
independent form factors 

+ - - + (S^v)}, (3.6) 



where 



('/if/2)"' " " (gi?2) 

From this amplitude the standard decay width follows 



r(Xco-77) = ^«L^^, (3.7) 

So the normalization of the amplitude (13.61) is correct. 

Analogously with the diffractive case (12.31) let us define the photon transverse momenta. 
Momentum conservation dictates us the following decompositions of photon momenta into 
longitudinal and transverse parts 

Qi = xipx + -jp 2 + <Jx,t, Q2 = -X2P2 ~ Pi + Q2,t, ql/ 2 ,t — -^1,2) , (3.8) 

where t\p = Qi2- Due to the gauge invariance we have similarly to (12.91) 

(Pi + p'lfV^M + p' 2 y = ip^V^. (3.9) 

In the relevant limit t\p <C xi, 2 s we get finally the following matrix element for pp — * pp Xc(0) 
reaction via 7*7*-fusion 

u rr „ _- A ^e 2 V^(0) [N7 F 1 (t 1 ) Fx(t 2 ) 3M 2 ( gi ,g 2 ,) + 2^2 - -higL 
^ \ 3 J M V vrM t x t 2 (M* - £ - ^ 

where (qi^.t) = — \/tit 2 (l — x\)(l — x 2 ) cos$, and $ is the relative angle between photons 
(or outgoing protons). The amplitude (13.101) is purely imaginary, so there is no interference 
with diffractive process describing by purely real amplitude (I2.5p . This amplitude will be 
used in the following to calculate differential cross section. 



B. Equivalent Photon Approximation 

In the equivalent photon approximation (EPA) the total cross section for pp — > pxdtyp 
can be written as a convolution of EPA flux factors and the 77 — > Xc(0) resonant cross 
section 

a = J d Zl dz 2 ^(^)^(z 2 )a(77 - xMYj ■ (3.11) 
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The elementary cross section can be written in terms of partial decay width as 

4tt 2 

*(ry - Xc(O)) « w T Xc{0) ^6(M - M R ) . (3.12) 

R 

Let us introduce two new variables: 

= (3.13) 



M = ^sziz 2 . 

Now the cross section can be written as 

da 2M dn . . dn , . 47r 2 ,./„,„,, 

5^5M = 1(^3 3£<*> S| r *W<* - Mn) . (3.14) 

Integrating over invariant mass of the two photons we get 

da _ 2M R dn dn An 2 , , 

dx~ F " s( Zl + z 2 ) ~dz~} Zl) dz~ 2 {Z2) Mf*>®-" ' (3 ' 15) 



where now 



z x = l -x F + ^x 2 F + AMl/s, 
1 1 



z 2 = --x F + -\/x 2 F + AM 2 R /s 



(3.16) 



This equation is suitable to calculate distribution of Xc(0) in the Feynman variable x F . The 
analytical flux factors of photons in protons of Drees and Zeppenfeld [Il| are taken. We 
have also tried: 

/M^(*)=/*ft (9?+ i, M , *?»), ("7) 

where t = -{q 2 + zm 2 N )/{l - z). The decay width from PDG [42] is r Xc(0) ^ 77 = 0.26 26 1(T 5 
GeV. 



IV. POMERON-POMERON FUSION 

Above we have shown how to calculate the diffractive Xc(0 + ) meson production mecha- 
nism in a QCD-inspired approach. Often in the literature in order to describe the high-energy 
processes one uses a phenomenological object known as pomeron. Often a vector nature is 
prescribed to such an object, i.e. it is assumed that it couples to the nucleons or (quarks) 
via 7^ matrices, i.e. similarly as photon. The corresponding mechanism for exclusive Xc 
meson production is sketched in FigJHJ 

The amplitude for our exclusive process may be written as 

«^i,ii)(Pi+P;) y C^ XcJ ( A; i'^)(P2+p / 2 )^( S2 ,t 2 ) . (4.1) 

In the equation above Ar(sx/2, ti/2) are so-called Regge propagators. They can be written 
as: 

A R (s 1/2 ,t 1/2 )=r c ^(^j 5 (4.2) 
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ho ho 
FIG. 8: Pomeron-pomeron fusion mechanism of x c -meson production. 

-Ff (£1/2) are form factors describing helicity-preserving coupling of the pomeron to the nu- 
cleon. We take them to be identical to the Dirac electromagnetic form factors of the proton. 
Such a choice is justified by the phenomenology of elastic proton-proton scattering (see 
e.g. 34|). The pomeron coupling to the charm quark and/or antiquark is reduced compared 



to the coupling to the nucleon (Cp) by the factor 3 (three quarks in the nucleon versus 
single charm quark/antiquark) and extra factor r c (heavy quark interaction is weaker than 
light quark interaction). This parameter can be extracted from the inelastic interaction of 
J/ip with the nucleons in nuclei lq . In the Tevatron energy r c ~ 0.2. The parameters 



of the pomeron exchanges: Cp and 5 are taken from the Donnachie-Landshoff fits to the 
proton-proton and proton-antiproton total cross sections |4(|. The vertex V w ~* Xc can be 
obtained from the vertex V 1 * 1 *^* by replacing charm quark charge by unity. 

The comparison of the results with pomeron-pomeron fusion with the diffractive QCD 
results will be given in the result section. 

V. CROSS SECTION AND PHASE SPACE 

The cross section for the 3-body reaction pp — > ppXc can be written as 

da„ PXc = ^\M\ 2 -d 3 PS. (5.1) 
The three-body phase space volume element reads 

l3 DC ^P'l d?p' 2 d 3 PftI ,„ <4rf, / I T~) \ /r n \ 

d PS = 2E[{2,f2E' 2 {2,f2E M (2,Y ' ^ 6 {Pl + » " * " " P ^ " ( 5 ' 2 ) 

At high energies and small momentum transfers the phase space volume element can be 
written as 

d 3 PS = -^-dhdhd^d&d® 5 (s(l - &)(1 - 6) - M 2 ) , (5.3) 

2°7T^ 

where £1, £2 are longitudinal momentum fractions carried by outgoing protons with respect to 
their parent protons and the relative angle between outgoing protons $ G (0, 2ir). Changing 
variables (£1,^2) — * (xf,M 2 ) one gets 

d 3 PS = JL-dhdh — dXF d$ . (5.4) 

2 8 ^ 4 V^4 + 4(M2 + |P Mit |2)/ s 
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It is more convenient for lower (but still high) energy to use variable xf- However, at 
very high energies the cross section becomes too much peaked at xf ~ due to the jacobian 

j « 1 _> f5 5) 

y/x% + AM 2 /s 2M v ' ; 

and the use of rapidity y instead of xp is recommended. The phase space element in this 
case has the following simple form 

d 3 PS = -r^—dhdhdyd® . (5.6) 

2 8 7T 4 S 

If Xp is used then 



i 1 / 9 4M2 ^ 

^ Kl -# + T T i' (5 - 7) 

In the other case when the meson rapidity is used then 

£ lj2 «l-^exp(±y). (5.8) 
Now the four-momentum transfers in both proton lines can be calculated as 

= _^ _ , ( , 9) 

Only if = 1, ti,2 = ~P\/2f ^he l a ^ter approximate relation was often used in earlier 
works on diffractive production of particles. However, in practice £i i2 7^ and the more 
exact equation must be used. The range of t% and ti is not unlimited as it is often assumed. 
One can read off from Eq. (15.91) a kinematical upper limit for t 12 which is 



\2 

*i,2 < , m 2 v . (5.10) 



1.2 



p 



In practice these phase space limits become active only for \xp\ > 0.2. The lower limits are 
energy dependent but are not active in practice. 



VI. RESULTS 

In Ref.jH estimates of the integrated cross sections were given. We wish to concentrate 
on differential distributions. Before we show our results we wish to discuss uncertainties 
related to the KMR approach. 



A. Uncertainties in the KMR approach 

In the KMR approach only one effective transverse momentum is taken explicitly in 
their skewed unintegrated distributions. In the KMR prescription it is the minimum of the 
transverse momenta of the two gluons connected to the same proton line. In order to see the 
uncertainties related to such a choice we shall present results obtained with the following, 
equally arbitrary, choices: 
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1) Ql,t = min (9o,t> (lit), Qh = min (#L> el,*) 

2 ) Qi, t = max(g£ tJ g 2 t ), Q| )t = max(gg t , 



9i,t 



3) — 9i,ti ^2,* 

5) = (qtt + 9?, t )/2, g|, 



9o,t > 



(9o,t + &)/2 



As an example in Figj9]we show results for the Feynman xf distribution. As can be seen 
from the figure there are rather large uncertainties related to the choice of the effective 
transverse momentum. 




-1 
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-2 
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Xc(0 + ) 
W = 1960 GeV 
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FIG. 9: Effect related to the choice of effective transverse momentum in the KMR UGDF on 
distribution in Feynman variable xf in the KMR approach. The solid lines are for the choice 1 
(upper) and choice 2 (lower) The dashed line is for choice 3, the dotted line for choice 4 and the 
dash-dotted line for the last possibility. The calculation was done for the Tevatron energy W = 
1960 GeV. 



The estimate by Khoze, Martin, Ryskin and Stirling in Ref . [16| was done for one selected 
value of the hard scale in the KMR UGDF /i 2 = M 2 ( )/4. In Fig{10]we show the dependence 
of the differential cross section on the value of the scale. There is a sizeable effect (about 
a factor of two) on the cross section (see also Table 1 with integrated cross sections). For 
xp ~ 0, the smaller /z 2 the bigger the cross section. 
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FIG. 10: Scale effect on distribution in Feynman x F in the KMR approach (KMR UGDF, KMR 
vertex). The solid line is for fj? = M 2 and the dashed line is for fj? = M 2 /4. The calculation was 
done for the Tevatron energy W = 1960 GeV. 



In order to demonstrate sensitivity to the nonperturbative region of small transverse 
momenta in FigJTT] we show results with the KMR UGDF cut off for small values of the 
effective gluon transverse momenta Q\ < Q 2 cut - One observes a rather strong dependence on 
the value of the cut-off parameter. The bigger the value of the cut-off parameter the smaller 
the cross section. The results for the integrated cross section are summarized in Table I. 
Even for large values of the cut-off parameter {Q 2 cut ~ 1 GeV 2 ) sizeable cross sections are 
obtained. 

The estimate in Ref. 16[ was done assuming that the gluons in the gg — * x c (0 + ) vertex 
are on mass shell, which is exact only in the infinite meson mass limit. In this paper (see 
subsection IIB) we take into account the effect of gluon virtualities. In FigJT2]we show the 
role of the off-shell effects for the KMR UGDF. The off-shell effect leads to a reduction of 
the cross section by a factor of 2 - 5. 

Finally we show the influence of the off-shell effects on azimuthal-angle correlation func- 
tion. The off-shell effect in the matrix element reduces the cross section but does not change 
its shape. The shape of the distribution requires an extra comment. The distribution in 
azimuthal angle is very different from (1 + cos(2$)) as expected for one-step fusion (e.g. 
photon-photon, pomeron-pomeron fusion) to be discussed in detail in the next subsection. 
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FIG. 11: The effect of the cut on low gluon transverse momenta on distribution in Feynman xf 
in the KMR approach (KMR UGDF, KMR vertex). The solid line is for Q 2 cut = 0.26 GeV 2 , the 
dashed line for Q 2 ut = 0.5 GeV 2 , the dotted line for Q 2 ut = 0.8 GeV 2 and dash-dotted line for Q 2 ut 
= 1.0 GeV 2 . The calculation was done for the Tevatron energy W = 1960 GeV. 



B. Our approach 

In FigJTHwe present our distributions of the Xc(0 + ) mesons in the Feynman variable xp. 
We present results obtained with different UGDFs. Characteristic for central diffractive 
production all distributions peak at Xp ~ 0. Although all UGDFs give a similar quality 
description of the low-a; HERA data for the F 2 structure function, they give quite different 
longitudinal momentum distributions of x c (0 + ) at the Tevatron energy W = 1960 GeV. The 
UGDFs which take into account saturation effects (GBW, KL) give much lower cross section 
than the BFKL UGDF (dash-dotted line). Similar as in Ref.[24| rather small values of x's 
in formula ( 12.5ft come into the game here. Therefore the process considered here would help, 
at least in principle, to constrain the poorly known UGDFs. 

In FigfTB] we show distributions in the square of the four-momentum transfers (t\ or i 2 ) 
in the nucleons lines. Because they are identical we shall denote them da/dt for brevity. 
The distributions shown in the figure are peaked at small values of ti or t 2 . Slightly different 
slopes are obtained with different UGDFs. This demonstrates how purely known are UGDFs 
at present. The measurement of such distributions requires measuring forward protons (or 
antiprotons). In the case no measurement of forward nucleons is possible one should study 
events with low multiplicity (only particles from the decay of Xc(0 + )). This can be either 
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FIG. 12: Off-shell effect on distribution in Feynman xp in the KMR approach. The dashed line 
is for on-shell case and solid line is for off-shell case. The calculation was done for the Tevatron 
energy W = 1960 GeV. 

pairs of pions or pairs of kaons or pairs of photons or J ftp + 7. Unfortunately branching 



ratios to these channels are rather low 42 



In FigJT6] we show azimuthal angular correlations of the outgoing nucleons. Measur- 
ing such distribution experimentally requires identification of both nucleons in very for- 
ward/backward directions. This is not possible with the present Tevatron apparatus. We 
hope such a measurement will be possible with the final LHC instrumentation. All distribu- 
tions shown in Fig{16]are peaked at $ ~ 180° i.e. for the back-to-back kinematics. The fact 
that the distributions are not simple functions (sin$, cos$) of the relative azimuthal angle 
between outgoing nucleons is due to the loop integral in Eq. fl2.5p which destroys the depen- 
dence one would obtained with single fusion of well denned objects (mesons or reggeons). 

Another type of proton-antiproton correlation is shown in FigJlTl We show results for 
KL and BFKL UGDFs as well as for photon-photon fusion. Finally we show the two- 
dimensional distribution in ti and ti for the "fusion" of two vector pomerons. The (£1,^2) 
distribution obtained in the photon-photon fusion mechanism differs qualitatively from the 
all other distributions. One can see a strong enhancement of the cross section when t\ — > 
or t2 — > 0. Also the shape of the two-dimensional spectrum obtained for the fusion of 
phenomenological vector pomerons differs from those obtained in the QCD-inspired KMR 
mechanism (please note minimum of the cross section when t\ — > and/or £2 — * 0). 

In order to demonstrate the role of the gluon transverse momenta in UGDFs in FigJTBI 
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FIG. 13: Off-shell effect on distribution on azimuthal correlation function in the KMR approach. 
The dashed line is for on-shell case and solid line is for off-shell case. The calculation was done for 
the Tevatron energy W = 1960 GeV. 

we show distributions in Xp for the Gaussian UGDF for different values of the smearing 
parameter cq. We wish to notice here that all such UGDFs correspond to identical integrated 
GDFs. The smaller value of Oq the larger the cross section. This demonstrates that the main 
contributions to the cross section come from the region of very small transverse momenta 
of the t-channel gluons (see FigH|). This is cleary the region where nonperturbative effects 
are dominant. 

Let us turn now to the second mechanism of the exclusive Xc(0 + ) production - the two- 
photon fusion mechanism sketched in FigfSJ In Fig{19] we show corresponding distribution 
in the Feynman variable Xp. For comparison we show also the EPA result (dotted line) 
discussed in subsection IIIB. 

For pedagogical purpose in FigfSO] we show also the result of the calculation with vector- 
like pomerons as described in section IV. Rather smaller cross section is obtained compared 
to the more QCD-inspired calculation with UGDFs. The cross section depends strongly on 
the value of the parameters of the form factor describing coupling of the phenomenological 
pomeron to nucleons. The value of B = 15-20 GeV -2 is preferred from the proton-proton 
or proton-antiproton elastic scattering phenomenology. Compared to the calculation with 
UGDFs presented in FigJTH the distribution around xp = is much broader here. The 
thinner distribution in the case of the /^-factorization approach is due to the x-dependence 
of the UGDFs entering the basic formula. The values of x's (x\, X2, etc.) in the off-diagonal 
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FIG. 14: Distribution in Feynman x F for different UGDFs (KL-dashed, GBW-dotted, BFKL- 
dash-dotted) for off-shell (left panel) and on-shell (right panel) matrix element. The calculation 
was done for the Tevatron energy W = 1960 GeV. 

UGDFs change quickly with the Feynman variable xp which makes the peak at = 
much thinner than for the phenomenological pomeron exchange discussed here. 

Let us return for a while to t-distributions. In Figj21] we show single distributions for 
photon-photon and Pomeron-Pomeron fusion mechanisms in t = t\ = £2 (relevant for one- 
nucleon tagged case) obtained by projecting the two-dimensional distributions shown in 
FigflTl The electromagnetic component peaks at very small values of t due to the photon 
propagators. 

The results for the diffractive mechanism depend strongly on the details of the calcu- 
lation. Here we show only one distribution of diffractive component for easy reference. 
Although the diffractive component is subjected to much stronger absorption effects than 
the electromagnetic one, it is clear that the diffractive component dominates. 

In Fig[22] we show differential distribution in relative azimuthal angle for 7*7* and 
Pomeron-Pomeron fusion mechanisms. One sees a typical (1 + cos(2$)) dependence charac- 
teristic for one-step exchanges. These distributions are very different from those shown in 
FigJTBl for pQCD diffraction where the underlying mechanism is more complicated due to a 
two-step nature of the process (see FigHj). 

Above we have calculated only bare distributions. Those are subjected to absorption 
effects. The absorption effects are usually included by multiplying the bare distributions by a 
soft gap survival probability. The gap survival probability was estimated to S 2 (Tevatron) « 
0.05 and S 2 (LHC) ~ 0.025 jig . Absorption leads therefore to a large reduction of the bare 
cross section. In principle, the absorption effects may modify also shapes of differential 
distributions. This requires a further detailed analysis which goes beyond the scope of the 
present paper. 

As discussed in this paper there are huge uncertainties in estimating the cross section for 
exclusive Xc(0 + ) production, much larger e.g. than for the exlusive J/ip production where 
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FIG. 15: Distribution in t = t\ = t<i for different UGDFs. The notation here is the same as in 
FigHH The solid line is for the Gaussian distributions with <7o = 1.0 GeV and the second choice 
of the scales in Eq. jSlgj) . 

the corresponding amplitude can be related to the amplitude of J /if) photoproduction in 
ep collisions j43|. Only a real experiment may shed more light on the dynamics of Xc(0 + ) 
production mechanism. Measuring centrally produced Xc(0 + ), one forward nucleon (proton 
or antiproton) and imposing a condition of rapidity gap in the second hemisphere should 
allow, at least in principle, such a measurement. 

In Table I we have combined the cross section integrated over the whole available phase 
space for exclusive production of the x c (0 + ) mesons. 

In Table 1 we have assumed W = 1960 GeV, S 2 = 0.1 and BR(x c (0+W J/ij)+i) = 0.01. 
The recently measured value of the branching ratio is 1.3 ± 0.11 % [42J. The last column 
shows possible contribution of the x c (0 + ) — > J /if) + 7 decay to the exclusive production of 
J/if) if the soft photon cannot be correctly identified. We get typically 0.1-5 nb, i.e. much 
less than the direct J/if> photoproduction (43|. It would be interesting to calculate the 
contribution of Xc(l + ) - * J/^Pl an d Xc(2 + ) — > J/^fj, where the corresponding branching 
ratios are order of magnitude larger. On the other hand, within the approximations made 
in Ref.0] the cross sections for diffractive production of y c (l + ) and Xc(2 + ) vanish. It would 
be interesting to go beyond the approximations of Ref . [9j . 



The exclusive production of Xc has been reported recently by the CDF collaboration 45 
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FIG. 16: Distribution in Feynman relative azimuthal angle for different UGDFs The notation here 
is the same as in Figll41 The solid line is for the Gaussian distributions with a$ = 1.0 GeV and 
the second choice of the scales in Eq. (j2.34p . 

with the upper limit for the cross section 

Oexc{vV ->P+ J/i> + 7 + P) < 49 ± 18(stat) ± 39(sys) pb , (6.1) 

within the CDF experimental cuts. These numbers cannot be directly compared with our 
results in Table 1 that do not include the experimental cuts. Such a comparison requires a 
Monte Carlo type analysis as was done in Ref . j46| . 



C. Energy dependence 

Up to now we have concentrated on Tevatron energy. We expect some results for exclusive 
Xc production in not too distant future. It would be also interesting to measure exclusive Xc 
production at different energies. The obvious choices are RHIC and LHC in the near future. 
In Figj23] we show distributions in x F for different UGDFs for RHIC, W = 200 GeV (left 
panel) and LHC, W = 14000 GeV (right panel). Compared to Tevatron the distributions in 
xp for RHIC are wider and distribution for LHC are more narrow, concentrated around xp 
= 0. While at RHIC energy different UGDFs give relatively similar results, at LHC energy 
there is a difference of several orders of magnitude between results obtained with different 
UGDFs. Therefore a measurement at LHC should clearly select the best distribution. 
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FIG. 17: Two-dimensional maps in h and t 2 for KL UGDF (top left) and BFKL UGDF (top right) 
as well as for two-photon fusion (bottom left) and two-pomeron fusion with dipole form factors 
(bottom right). 

In Table II we have collected total cross sections for selected UGDFs at RHIC, Tevatron 
and LHC. Comparing results for the three different energies we see that different UGDFs 
predict completely different energy dependence. While BFKL predicts a strong growth of 
the cross section with the collision energy, the saturation models (KL, GBW) predict much 
slower rise with the collison energy. Therefore measuring exclusive Xc(0 + ) production at 
three different energies would be useful to pin down underlying dynamics. 

VII. DISCUSSION AND CONCLUSIONS 

In the present paper we have discussed in detail the exclusive production of the x c (0 + ) me- 
son in proton-antiproton and proton-proton collisions. We have considered both diffractive 
and purely electromagnetic mechanisms. Many differential distributions have been discussed 
for the first time in the literature. 

The diffractive component was calculated in the Khoze-Martin-Ryskin approach. Com- 
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TABLE I: Integrated cross section (in nb) for exclusive Xc(0 + ) production for the Tevatron energy 
W = 1960 GeV in various possible theoretical prescriptions. The number for the Gaussian UGDF 
refers to the number in Eq, (|2.34p . 



UGDF, details 



Otot 



S 2 a tot 



BR S 2 a tot 



1. KMR, 
vertex, 

2. KMR, 
vertex, 

3. KMR, 
vertex, 

4. KMR, 
vertex, 

5. KMR, 
vertex, 

6. KMR, 
vertex, 



Q 2 cut = 0.26 GeV 2 , on-shell 
min prescription for Q 2 

Q 



cut 



0.26 GeV 2 , on-shell 
max prescription for Q 2 
Q 2 cut = 0.26 GeV 2 , off-shell 
min prescription for Q 2 
Ql ut = 0.5 GeV 2 , on-shell 
min prescription for Q 2 
Q 2 cut = 0.8 GeV 2 , on-shell 
min prescription for Q 2 
Q 2 cut = 1.0 GeV 2 , on-shell 
min prescription for Q 2 



0.1357(+4) 
0.9628(+4) 
0.3720(+4) 
0.8227(+3) 
0.4124(+3) 
0.2745(+3) 



0.1357(+3) 



0.9628(+3) 



0.3720(+3) 



0.8227(+2) 



0.4124(+2; 



0.2745(+2) 



0.1357(+1) 



0.9628(+l) 



0.3720(+l) 



0.8227(+0) 



0.4124(+0) 



0.2745(+0) 



4. KL, on-shell vertex 

5. KL, off-shell vertex 



0.1180(+3) 
0.5231(+2) 



0.1180(+2) 
0.5231(+1) 



0.1180(+0) 
0.5231(-1) 



6. GBW, on-shell vertex 

7. GBW, off-shell vertex 



0.8514(+2) 
0.1590(+3) 



0.8514(+1) 
0.1590(+2) 



0.8514(-1) 
0.1590(+0) 



8. BFKL, on-shell vertex 

9. BFKL, off-shell vertex 



0.2603(+4) 
0.1125(+4) 



0.2603(+3) 
0.1125(+3) 



0.2603(+l) 
0.1125(+1) 



10. Gauss, o"o = 0.5 GeV, 
off-shell vertex, scales (2) 

11. Gauss, o"o = 1.0 GeV, 
off-shell vertex, scales (2) 



0.2141(+2) 
0.1811(+1) 



0.2141(+L 



0.1811(+0) 



0.2141(-1) 



0.1811(-2) 



TABLE II: Integrated cross section otot (hi nb) for exclusive Xc(0 + ) production at different energies. 
The number for the Gaussian UGDF refers to the number in Eq. (|2.34j) . 



UGDF 


RHIC 


Tevatron 


LHC 


Kl 


0.6430(+l) 


0.5231(+2) 


0.1090(+3) 


GBW 


0.2830(+l) 


0.1590(+3) 


0.1413(+3) 


BFK1 


0.6140(+1) 


0.1125(+4) 


0.6306(+5) 


Gauss, (To = 1.0 GeV, 
scales (2) 


0.7126(-1) 


0.1811(+1) 


0.1428(+2) 
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FIG. 18: Distribution in xp of the Xc(0 + ) meson for the Gaussian UGDF. In the left panel different 
values of the parameter gq = 0.5, 1, 2 GeV and /i§ = fi 2 = M 2 . In the right panel do = 1 GeV 
and we have chosen different scales in Eq. (|2.34p : 1: solid, 2: dashed and 3: dotted. 



pared to the original KMR calculation we have taken into account the off-shellness of the 
gluons in the corresponding diagram. The corresponding matrix element was calculated. 
We find that the inclusion of the gluon virtualities reduces the cross section by a factor of 2 
- 5, depending on kinematical region and on UGDFs. 

We have discussed uncertainties in the KMR approach related to the treatment of the 
nonperturbative region as well as related to the choice of the scale in their skewed uninte- 
grated distributions. This gives an uncertainty of a factor 2-3. 

We get similar integrated cross section as in Ref. 16| if we take the on-shell approximation 
for the vertex, make the same choice of scales, use the same integrated gluon distributions, 
etc. Summarizing, we find rather large uncertainties in the approach. 

Many other UGDFs from the literature were also used to calculate distributions in the 
Feynman variable Xp, in squares of the four-momentum transfers (ti and t 2 ) in the nucleon 
lines and in relative azimuthal angle between outgoing protons. Also correlations in t\ and 
£2 have been analyzed. The results depend strongly on the choice of UGDF. This is related 
to a particular sensitivity of the cross section for the reaction under consideration to the 
nonperturbative region of very small gluon transverse momenta. Therfore a measurement of 
the differential distributions would be very helpful to test the unintegrated distributions in 
this region. At RHIC one could measure the n + ir~ and/or K + K~ decay channels. Probably 
ALICE could use a similar method. At Tevatron rather 77 or 'yJ/ip channels seem preferable. 

In the present paper we have calculated only bare distributions. Those are subjected 
to absorption effects. The absorption effects are usually included by multiplying the bare 
distributions by a soft gap survival probability. The survival probability was estimate as 
S 2 (Tevatron) « 0.05 and S 2 (LHC) ~ 0.025 Absorption leads therefore to a large 

reduction of the bare cross section. In principle the absorption effects may modify also 
shapes of differential distributions |43j. This requires a further detailed analysis which goes, 
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FIG. 19: Distribution in xp for the two-photon fusion. The solid line represents result obtained 
with Eq. (|3.10p . The dotted and dashed dotted lines correspond to the EPA approach with different 
flux factors as described in the text. 

however, beyond the scope of the present paper. 

For completeness we have calculated also the cross section in a more phenomenological 
approach with Regge-type pomeron-pomeron fusion. Cross sections of the same order of 
magnitude as for the QCD approach are obtained. However, the azimuthal angle correlation 
functions for both approaches are very different. 

We have also calculated differential distributions for the photon-photon fusion. This con- 
tribution turned out to be rather small (a fraction of nb). The differential distributions have 
quite different shapes compared to the diffractive component. Compared to the diffractive 
component it is peaked at extremely small values of t\ and/or ti- Also the azimuthal angle 
correlations pattern is very different: cos$ for the 77 fusion and a more complicated shape 
for the diffractive component due to the loop integration in the formula for the amplitude. 

Recently, there is an interest at Tevatron to study exclusive production of J_hp meson. 
This is because of a possibility to study J/ip pho toproduction at large energies [43|] and due 
to a potential search for odderon exchange [44|. Because of its decay channel Xc(0 + ) — > 
•yJ/ip and difficulties in identifying soft photons the Xc(0 + ) decays may contribute to the 
"exclusive" production of J/\l/ making the measurement of J/^f photoproduction and/or the 
discovery of the odderon exchange difficult. A comparison of theoretical cross sections for the 
three mentioned reactions would be therefore very useful in this context. Our estimates for 
Tevatron energies, including soft survival probability, give (20 - 200) nb times BR(xc{^ + ) — ► 
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FIG. 20: Distribution in xf for the two-pomeron fusion for different values of the slope parameter 
B = 5, 10, 15, 20 GeV -2 (from top to bottom) of the exponential form factor. The dashed line is 
obtained with the dipole (electromagnetic) form factor. 

j J /if}) = 0.013 which gives the corresponding cross section rather less than 1 nb compared 
to about 15 nb for the J /if) photoproduction (43| . 

We have made also calculation of da/dxF distributions and total cross sections for RHIC 
and LHC Comparing these results and those for the Tevatron leads to the conclusion that 
only measurements for different energies may help to disentangle underlying QCD dynamics. 
In principle, such measurements could even explore the onset of QCD saturation which is 
not so easy to be discovered in inclusive reactions. 
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